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I. INTRODUCTION 



The mechanism of Electro- Weak Symmetry Breaking (EWSB) will be one of the primary 
topics to be investigated at the LHC, both through direct searches for the Higgs boson and 
careful analysis of boson boson scattering processes. Detailed reviews and extensive bibli- 
ographies can be found in Refs. . The nature of the interaction between longitudinally 
polarized vector bosons and the Higgs mass, or possibly the absence of the Higgs particle, 
are strongly related: if a relatively light Higgs exists then the Vl's are weakly coupled, while 
they are strongly interacting if the Higgs mass is large or the Higgs is nonexistent 

It should be noted that the Goldstone theorem and the Higgs mechanism do not require 
the existence of elementary scalars. It is conceiveable and widely discussed in the literature 
that bound states are responsible for EWSB. 

At the LHC no beam of on shell EW vector bosons will be available. Incoming quarks will 
emit spacelike virtual bosons which will then scatter among themselves and finally decay. 
These processes have been scrutinized since a long time in order to uncover the details of the 
EWSB mechanism in this realistic setting fly , |7| . Naively one expects that at large boson 
boson invariant masses the diagrams containing vector boson fusion subdiagrams should 
dominate the total cross section while the offshellness of the bosons initiating the scattering 
process should become less and less relevant. Together with the keen interest for Higgs 



production in vector boson fusion, this has lead to the development of the Equivalent Vector 
Boson Approximation (EVBA) ja, 2I, 10]. The EVBA provides a particularly simple and 
appealing framework in which the cross section for the full process is approximated by the 
convolution of the cross section for the scattering of on shell vector bosons times appropriate 
distribution functions which can be interpreted as the probability of the initial state quarks 
to emit the EW bosons which then interact. 

This approach relies on the neglect of all diagrams which do not include boson boson scat- 
tering subdiagrams and on a suitable on-shell projection for the scattering set of diagrams. 
Since the approximate boson boson interaction is expressed in terms of on shell particles it 
is straitforward to separate the different boson polarizations. It is well known that the set 
of scattering diagrams is not separately gauge invariant while both the on shell amplitude 
and the distribution functions which appear in the EVBA are gauge independent. 

In Uj it has been shown that when vector bosons are allowed to be off mass shell in boson 
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bosoa scattering, the amplitude grows faster with energy compared with the amplitude for 
on shell vectors. Subsequently, it has been pointed out in [12| that the problem of bad high 
energy behaviour of WW scattering diagrams can be avoided by the use of the Axial gauge. 



Recently WZ production at hadron collider 



13] has been analyzed in Axial gauge with very 



encouraging results. It should be noted that the results in [12| have been obtained under the 
assumption that the transverse momenta of the produced W's are of the order of the Higgs 
mass and that each are much larger than the W mass. It is therefore not obvious to what 
degree the conclusions of Ref. 12) can be applied in the LHC environment, particularly for 
light Higgs masses as preferred by global SM fits. 

The EVBA has some undesirable features: a number of unphysical cuts need to be 
introduced to tame the singularities generated by the onshell projection, which are absent 
from the exact amplitude. In the literature, a number of comparisons of exact and EVBA 
calculations have appeared with conflicting results 

For these reasons in the present paper we have critically examined the role of gauge 
invariance in VV-fusion processes and the reliability of the EVBA in describing them. We 
would like to determine regions in phase space, at least in a suitable gauge, which are 
dominated by the scattering set of diagrams. 

If this was the case then it would be possible employ this set of diagrams to define a 
"signal", that is a pseudovariable which could be used to compare the results from the 
different collaborations. The signal is not necessarily directly observable but it should be 
possible to relate it via Monte Carlo to measurable quantities. If such a definition is to be 
useful it must correspond as closely as possible to the process which needs to be studied 
and the Monte Carlo corrections must be small. If instead, as we are led to believe by the 
results shown in the following, the VV scattering diagrams do not constitute the dominant 
contribution in any gauge or phase space region, there is no substitute to the complete am- 
plitude for studying boson fusion processes at the LHC. We cannot claim to have examined 
all possible gauge scheme, however we have studied the most commonly used £-type gauges 
and the implementation of the Axial gauge proposed in 12|. The results obtained with the 
full amplitude show that the mechanism of EWSB can indeed be investigated even without 
separating out the scattering set of diagrams. 
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II. VV -SCATTERING AND GAUGE INVARIANCE 



In order to study the implications of gauge invariance in VV-fusion processes we have 
concentrated on the specific process us — > cdW + W~ in proton proton collisions at the 
LHC. The corresponding Feynman diagrams can be classified in two different sets (Fig. [T]), 
the boson boson fusion one (Fig. [T^l) and all the rest (Fig. [lb) in which at least one final 
W is emitted by a fermion line. The two sets are not separately gauge invariant. We 
have started our analysis from the contribution of the boson boson fusion diagrams and 
their interference with the remaining ones. We have considered the Unitary, Feynman and 
Landau gauge. Among these, we find that the Feynman gauge is the one which minimizes 
the cancellations between set (a) of Fig. [T] and the non scattering diagrams. We have also 
petfood the calculation in the Axia. gauge = withiu the scheme proposed iu 

Ref.[12|. In Appendix. A we have collected the main corresponding Feynman rules. We 
have tried different n M gauge vectors and we find that the choice = (1, 0, 0, 1) (The 
incoming protons propagate along the z axis) is the best one. For brevity we will refer to this 
framework as Axial gauge in the following. We present results for the Unitary, Feynman and 
Axial gauge, showing the contribution of all diagrams (all), of the WW scattering diagrams 
(WW) together with their ratio WW/ all. When this ratio is significantly greater than one 
the contribution of non scattering diagrams is of the same order as the contribution of the 
WW scattering ones and important cancellations take place between the two sets. 

In Tab. [I] the total cross sections and their ratios are computed in the limit of infinite 
Higgs mass. This limit will also be referred to in the following as noHiggs case. We find that 
the Axial gauge is the one in which the interferences are least severe with a ratio of about 
2, while the ratios for the Unitary and Feynman gauge are 358 and 13 respectively. The 
inclusion of a light Higgs (M/j=200 GeV) does not improve matters: on the contrary the 
ratios become even larger as shown in Tab. [Ill for M^=200 GeV, doubling for the Unitary 
and Feynman gauge 

The comparison of total cross sections is only a preliminary step. It provides the general 
behavior but it cannot give informations on the different regions of phase space. For this 
reason we have evaluated the distribution of different kinematical variables with the goal of 
finding regions where the interferences are not important. If such regions exist it will be 
possible to define a set of cuts which allow the extraction of the WW scattering amplitude. 
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(a) 




(b) 



FIG. 1: Main diagram topologies for the process us — > cdW + W 





cr(pb) All diagrams 


cr(p6) WW diagrams 


ratio WW /all 


Unitary gauge 


1.86 10~ 2 


6.67 


358 


Feynman gauge 


1.86 1(T 2 


0.245 


13 


Axial gauge 


1.86 10~ 2 


3.71 1(T 2 


2 



TABLE I: WW diagrams and complete set of diagrams cross sections and their ratios computed in different 
gauges without Higgs contribution. We have used the CTEQ5 Pdf set with scale M\y 

Fig. [2] shows the distributions of the diboson invariant mass for the complete set of diagrams 
and for the WW diagrams only together with the ratio (da(WW) / dM) / (da(All) / dM) . We 
have chosen the vector boson pair invariant mass distribution as a prototype but the same 
conclusions are reached with all other variables. The results of Fig. [2] have been obtained for 





cr(pb) All diagrams 


cr(pb) WW diagrams 


ratio WW 1 all 


Unitary gauge 


8.50 10~ 3 


6.5 


765 


Feynman gauge 


8.50 10~ 3 


0.221 


26 


Axial gauge 


8.50 10~ 3 


2.0 10~ 2 


2.3 



TABLE II: WW diagrams and complete set of diagrams cross sections and their ratios computed in different 
gauges with a M h =200 GeV Higgs and M{WW) > 300 GeV . 
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a very large Higgs mass but the general behaviour is not modified by the inclusion of a light 
Higgs. The distribution obtained with the full set of diagrams and with the WW fusion set 
only are quite different in Unitary and Feynman gauge. The ratio is also large over the whole 
interval and especially in the region of high invariant mass which is the most important one 
for EWSB studies. Again the Axial gauge gives the best result, with distributions which 
have the same general shape in the two cases. The ratio in this gauge remains however 
greater than 2, apart from a very small region at low invariant mass. 



WW invariant mass distribution t 




M„ w (GeV) Www (GOV) 



FIG. 2: Distribution of da/dMww for the process PP — » us —* cdW + W~ for All diagrams, WW 
diagrams and their ratio in Unitary, Feynman and Axial gauge in the infinite Higgs mass limit. 
The Unitary gauge data in the left hand plot have been divided by 20 for better presentation. 



We have completed our study by analyzing bi-dimensional distributions of several pairs 
of kinematical variables. These double distributions allow us to analyze in particular the 
situation where the two incoming bosons have small virtualities (tip — ► 0). In this region the 
WW scattering diagrams are expected to dominate. The same behaviour found previously 
is observed: the effect of the interferences is much more relevant in Unitary gauge with 
respect to the others. In general the WW fusion subset has a different behavior compared 
with the complete calculation. In Fig. [3] we report dajdtidti in the Axial gauge where 
tip = tj—(pu,s — Pd,c) 2 are the square root of the absolute value of the invariant masses of 
the incoming off shell W's. The corresponding ratio distributions (plots (c) (d)) show that 
even in the limited region < tip < 250 GeV, the WW contribution alone does not provide 
a realistic description of the complete set of diagrams. Even in Axial gauge significant 
cancellations take place. The ratio is highly asymmetric in the ti — t 2 plane, reflecting its 
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sensitivity to the choice of the Axial gauge axis (Recall that in this case the gauge vector 
is along the z-axis) and WW scattering diagrams alone provide only a very rough estimate 
of the cross section. It is however possible to find regions, typically for t li2 > 200 GeV 
where the cross section is greatly reduced, where the ratio of the two results is around one. 
For comparison we also show the distributions obtained from the WW fusion subset in 
the Unitary (e) and Feynman (f) gauge whose shape, not to mention the normalization, is 
completely different from the results obtained from the full amplitude. 

From this sample of results one can conclude that the WW scattering diagrams do not 
constitute the dominant contribution in any phase space region for the gauges we have 
examined. For the Axial gauge, which in various cases shows ratios of order 2, we remark 
that these have been obtained only with a particular time like gauge vector = (1, 0, 0, 1). 
For any other space-like n^, the interferences are so important that even the numerical 
integration becomes difficult for the WW diagrams set whereas no problem occurs for the 
full set of diagrams. So we did not succeed in finding a gauge vector of the type suggested 
in 12| for which the cancellations are negligible. Even with our best choice of the 
distributions of the different kinematical variables show the presence of large interferences 
between the two subset of diagrams over most of phase space, indicating that the WW 
diagrams are not dominating. As a consequence a question mark is put on the possibility to 
isolate the WW scattering contribution by restricting the calculations to the corresponding 
diagrams. The reliability of approximation methods based on such an approach becomes 
then suspicious. For this reason we have completed our analysis by studying the Effective 
Vector Boson Approximation applied to our prototype process. 
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FIG. 3: Double differential distributions of da/dt\dt2 (pb/GeV 2 ) with t± t 2 = y—(p u ,s — Pd.c) 2 for 
the process PP — > us — > dcW + W~ , using All diagrams (a) or only the WW fusion subset in Axial 
gauge (b). The two plots in the central row represent the ratio WW /all in the Axial gauge as a 
function of t\,t2 for the total region (c) and for small t's (d). For comparison we also show the 
distributions of the WW fusion subset in the Unitary (e) and Feynman (f) gauge. All invariants 
are in GeV. 8 



III. THE EFFECTIVE VECTOR BOSON APPROXIMATION 



The Effective Photon Approximation known also as Weizsacker- Williams approximation 



IS 



171 ] has proved to be a useful tool in the study of photon-photon processes at e + e~ 
colliders. Encouraged by this success the approach has been extended to processes involving 
massive vector bosons 0, [h]] under the name of Effective Vector Boson Approximation. 
The EVBA has been first applied at hadron colliders in connection with Higgs production 



pp — > H + X 18] and subsequently to vector boson processes of the type pp — > (H — > 
V3V4) + X — > V3V4 + X 0, Q] in order to obtain EVBA predictions for the production 
of a vector boson pair not necessarily near the Higgs resonance and to study the strongly 
interacting scenario. 

In analogy to the QED case, the application of EVBA to these processes consists in: 

1. Restricting the computation to the vector boson scattering diagrams neglecting dia- 
grams of bremsstrahlung type. 



2. Projecting on-shell the momenta of the vector bosons which take part to the scattering: 
f 2 

l Vl, 2 



q\ 2 = My . Here it is important to notice that contrary to the 77 processes where the 



photon momentum can reach the on shell value qf 2 = 0, for the massive vector bosons 
the onshell point qf 2 = My 2 is outside the accessible phase space region q\ 2 < 0. 

3. Approximating the total cross section of the process fif 2 —* /3/4V3V4 as the con- 
volution of the vector boson luminosities Cp d f Po i 2 (x) with the on shell vector boson 
scattering cross section: 

<t(/i/ 2 -> /3/4V3V4) = / £ E C^p^a^V^V^xs^dx (1) 



V\,V% PohPoh 



Here x = M(V\V2) /s qq , while MiVilty is the vector boson pair invariant mass and s qq 
is the partonic center of mass energy. 

It is clear that this approximation provides a simplification from the computational point 
of view and can exploit the properties of on shell boson boson scattering. In first applications 
of EVBA to vector boson scattering further approximations have been adopted. The zero 
angle scattering approximation has been used for the process V1V2 — > V3V4 so the transverse 
momentum of the incoming bosons has been neglected. The zero mass limit for the vector 
boson mass My — > has also been considered in the computation of the luminosity. In this 



early application of EVBA only the contribution from longitudinal modes was considered 
while the contribution from transverse states was neglected. The EVBA results depended 
strongly on the details of the approximations made. The EVBA generally overestimated 



the complete perturbative calculation, in some cases by a factor 3 [20|. Progressively, more 
refined and rigorous formulations of EVBA have been proposed, avoiding as much as possible 
the mentioned approximations, until a formulation where no kinematical approximations 



are taken and a 
into account 



1 vector boson polarization states as well as their interferences are taken 



2l| . Inspired by the strategy of (2lj, we have further improved the EVBA. 
We have used an approach where not only all kinematic approximations are avoided but 
also the luminosity computation is not needed. This allows in particular to keep all final 
particle properties (momenta, angles...) as they would be in an exact calculation contrary 
to the traditional approach where a pre-integration is performed to obtain the vector boson 
luminosities. 

IV. EVBA APPLICATION TO THE PP -» us -> dcW + W~ PROCESS 



In 2l| a precise formulation of EVBA has been developed. It is based on a factorization 



technique for analyzing Feynman diagrams which leads to exact probability distribution 



functions for the vector bosons [22[. This improved formulation does not invoke any kine- 
matical approximation such as those mentioned in the previous paragraph. The only ap- 
proximation concerns the on shell continuation of the vector boson scattering cross section. 
Using this factorization technique and the relation between the polarization vectors and the 
vector boson propagator in Unitary gauge, the matrix element for any process of the kind 
fif 2 — > ^3/4 + Y, where Y is produced by vector boson fusion, can be written as: 

MS t x x M K „) (2) 

m,n=-l 1l ~ 1V1 Vi <?2 - M V 2 

qip are the momenta of the initial vector bosons, 6j(m) are their polarization vectors 
corresponding to the different helicity states m = 0, ±1. We have used the same expressions, 
conventions and frame to define them as given in [2l{|. They are normalized according to : 

e J (m)-e*(m')=5 m , m i-l) m (3) 
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and satisfy the completeness relation: 



E WH-f + S (4) 

m=- 1,0,1 m Vj 

ji t 2 are the quark currents and M(m, n) the off shell scattering amplitude of the vector 
boson subprocess V1V2 — > V3V4: 



M(m,n) = £ <£(roK(n)T^e$°ef (5) 

Usually an integration is performed over all the integration variables which are not con- 
cerned with the vector boson scattering subprocess as a first step to compute the vector 
boson luminosities. Denoting these variables by {0}, one can write the total cross section 
expression as: 

otoi = J g(q 2 , ql, x, (f))d(f) a^ f (W 2 , q 2 , q 2 ^)dq\dq\dx (6) 

g(qf, q%, x, 0) represents all terms which are independent of the vector boson scattering 
subprocess. Here x = — , where W 2 = (q\ + 02) 2 is the diboson invariant mass squared. 

Sqq 

a V J } (W 2 ,qtql)~ j £ M(m,n)M*(m' ,n')dp v (7) 

m, n, m'n' 



is the off shell VV — > VV cross section and dpy the final state vector 



element. The next step is the extrapolation to on shell masses. In 2l| it is achieved by 



Doson phase space 



simple proportionality factors between the off shell and on shell cross sections: 

° M(W 2 , ql q 2 2 ) = f vol (W\ ql q\)^ ol {W 2 \ M 2 ^ M 2 2 ) (8) 
the subscript pol refers to the different vector boson polarization states, and M^, My the 



masses of the vector bosons initiating the scattering process. While we refer to [21y for 
the details, we reproduce here the form factors f po i expression according to the different 
polarization configurations: 



M 2 M 2 

Ar = l; t = (9) 

— Qi —q2 

M 2 M 2 

hr = A) ; hL = A) (10) 

— yi — <?2 
II 



TLTL 



M Vl My 2 



(ii) 



Finally one can write the cross section expression in the EVBA approximation as: 



a 



EVBA 
tot 



M, ?2> x > 4>)d(j)dq{dql ^ fvoi^dx 

PohPol 2 



(12) 



By comparison with eq.((T]) the luminosity can be expressed as: 



^Zpo1 2 ( x ) = / g(qf,qix,(i>)fpoid(j)dqldq% 



(13) 



In our implementation we have used the same assumptions but we have performed the 
onshell extrapolation at the matrix element level. This allows to keep all the terms 
M(m,n)M*(m f ,n') (m ^ m',n ^ n') in the total amplitude square expression. The off- 
shell vector boson scattering matrix element M(m, n) in ([2]) and ([5]) is expressed in terms 
of the corresponding on shell matrix elements and the polarization factors (Q as 



M(m, n) = JfZM on (m } n) (m, n = L,T). 



(14) 



Moreover we have not employed any luminosity function but we have used the diagram- 
matical expression of the fermion lines, evaluating them for each kinematical configuration 
and polarization. 



Invariant mass distribution of the W pair | 




Invariant mass distribution of the W pair 




200 300 400 500 600 700 800 900 1000 
Mww (<»ev) 



200 300 400 500 600 700 800 900 1000 
Mww (QeV) 



FIG. 4: WW invariant mass distribution M(WW) for the process us — > dcW + W with EVBA 
(black solid curve) and with exact complete computation (red dashed curve) for the Infinite Higgs 
mass case (left) and M h =250 GeV (right) at the LHC. 
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V. COMPARISON BETWEEN EVBA AND EXACT RESULTS 



We have computed the total cross section and the distribution of the vector boson pair 
invariant mass M(WW) with an exact calculation and in EVBA for the process PP — > 
us — > cdW + W~ at the LHC. Cuts are necessary in EVBA to avoid the photon t channel 
propagator pole and the form factor (Q (qf — > 0, q\ — > 0) singularities. We have restricted the 
CM scattering angle of the W's and the laboratory frame polar angle of the c and d quarks in 
order to deal with photon and form factor singularities respectively. For comparison purpose, 
the same cuts have also been applied in the exact calculation even though no singularity 
appears in this case. Unless otherwise noted, we have used a cut of 10 degrees for all three 
angles. 

For the total cross section in the noHiggs case we obtain 0.63 10~ 2 pb with an exact 
calculation and 1.36 10~ 2 pb in EVBA. For 714=250 GeV we obtain 1.32 10~ 2 pb and 1.64 
10 -2 pb respectively. The corresponding distributions of the WW invariant mass is shown 
fig. HI The ratio of the two results is . The EVBA overestimates the exact calculation by 
a factor of about two, which is almost insensitive to M(WW), with the exception of the 
Higgs peak region at M^= 250 GeV where the exact result is larger than the EVBA one. 

We have checked our implementation of the EVBA against the one used in PYTHIA for the 
scattering of longitudinally polarized W's. The two results are only in qualitative agreement 
as expected both for the total cross section and for distributions. 



M h 


EVBA (pb) 


EXACT (pb) 


Ratio 


oo 


3.90 10" 2 


1.78 10~ 2 


2.17 


130 GeV 


3.94 10~ 2 


1.71 10~ 2 


2.3 


250 GeV 


4.61 10~ 2 


4.09 10~ 2 


1.12 


500 GeV 


4.42 10" 2 


2.5 10~ 2 


1.77 



TABLE III: Total cross sections computed with EVBA and exact computation and their ratio for the 
process us — > cdW + W~ at fixed CM energy y/s — 1 TeV. 



In previous works 7j the same partonic process us — > cdW + W~ has been used for EVBA 
vs exact comparisons at fixed energy. For this reason we have tried to reproduce the results of 
[zj], using as closely as possible the same cuts and parameters. On one side this represented 
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@cut 


EVBA (pb) 


EXACT (pb) 


Ratio 


10° 


4.42 1(T 2 


2.5 10~ 2 


1.77 


30° 


1.33 10~ 2 


2.06 10~ 2 


0.64 


60° 


6.06 10~ 3 


1.28 10~ 2 


0.47 



TABLE IV: Total cross section in EVBA and exact computation and their ratio for different angular cuts. 
The CM energy is ^/s = 1 TeV and the Higgs mass M h =500 GeV. 

a test of our implementation of EVBA and on the other side it allowed us to determine 
whether the EVBA results are improved by this new and more rigorous formulation. We 
have computed the total cross section at a fixed center of mass energy of 1 TeV in the 
limit of infinite Higgs mass and with Higgs masses M^=130, 250, 500 GeV. In Tab. IHII the 
cross sections and corresponding ratio are reported. We have also compared a number of 
distributions obtained with the two versions of EVBA and with the exact calculation. Our 
more sophisticated implementation doe not appear to improve substantially the agreement 
with the exact results. 

Finally, we have analyzed the sensitivity of the results to the angular cut in the M^=500 
GeV Higgs case. The total cross section for 9 cut = 10°, 30°, 60° are presented in Tab. IIVI 
which shows that the EVBA is more sensitive to the angular cut than the exact computation. 
The corresponding M(WW) distribution is shown in fig [H We see that the relationship 
between the exact and EVBA results depends quite appreciably on the angular cut. The 
EVBA overestimates the exact result at 9 cut = 10° by about a factor of two outside the 
Higgs peak, while the two distributions are in fair agreement at the resonance. However 
the EVBA underestimates the correct result at 8 cut = 60° over the whole mass range. The 
difference decreases from about a factor of two at small invariant masses to roughly 20% 
at masses larger than the Higgs mass. Therefore, while it appears quite possibile to find a 
set of cuts, at fixed energy and Higgs mass, for which the EVBA approximation reproduces 
well the exact result for the total cross section, in general it is extremely difficult to extract 
from the EVBA more than a very rough estimate of the actual behaviour of the Standard 
Model predictions for boson boson scattering. 
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Invariant mass distribution of the W pair 




200 300 400 500 600 700 800 900 

Mww (GeV) 

FIG. 5: The WW invariant mass distribution with different #-cut with EVBA (solid curves) and 
with the exact calculation (dashed curves). From top to bottom 9 C = 10° (black), 30° (red), 60° 
(blue). The CM energy is y/s = 1 TeV and M h =500 GeV 
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W+W- Invariant mass distribution 




1000 1500 2000 2500 

Mww (GeV) 



FIG. 6: The WW invariant mass distribution in PP — > us — > cdW + W~ at the LHC for the infinite 
Higgs mass case (solid curves) and for Mh=200 GeV (dashed curve). The two intermediate (red) 
curves have been obtained imposing the set of cuts described in the text. The two lowest (blue) 
curves refer to the process PP — » us — > cdfj,~V fJj e + i> e ; in this case further acceptance cuts have 
been imposed on the charged leptons: E[ > 20 GeV, px t > 10 GeV , \rji\ < 3. 



VI. THE LARGE INVARIANT MASS REGION 

The results presented in Sects. HI} IVllead to the conclusions that the boson boson subam- 
plitude cannot reliably be extracted from the full qq — > qqVV amplitude. However the full 
amplitude is sensitive to the details of the EWSB mechanism. If no light Higgs is present in 
the SM spectrum, some hitherto unknown mechanism must intervene to enforce unitarity of 
the S-Matrix which embodies the conservation of total probability. The infinite Higgs mass 
limit violates perturbative unitarity in onshell boson boson scattering but on the other hand 
can be computed exactly, while the many available models for unitarizing the theory deal 
exclusively with on shell bosons and can only approximately be incorporated in a decription 
of qq — > qqVV processes or in a six final state fermion framework, which has recently become 
available for the LHC 



23 



24]. 



At the LHC the linear rise of the cross section at large boson boson invariant masses 
squared entailed by the leading behaviour of boson boson scattering in the SM with a very 
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E (quarks) > 20 GeV 
p T (quarks, W)> 10 GeV 
2 < |r?(quark)| < 6.5 
MW)| < 3 

TABLE V: Selection cuts applied in Fig. [6l 



large Higgs mass will be overcome by the decrease of the parton luminosities at large x and 
will be particularly challenging to detect. In the absence of a more reliable theory we have 
adopted the noHiggs model as a poor man substitute. 

In Fig. [6] we show the large mass tail of the boson pair invariant mass distribution for 
the noHiggs case and for a Higgs mass of 200 GeV. In the absence of cuts the two results 
differ by about 20% over the full range. With appropriate cuts the difference between the 
two cases can be significantly increased. Applying the selection cuts in Tab. |V] we obtain 
the two intermediate, red curves in Fig. [61 The two lowest (blue) curves refer to the full 
process PP — > us — > cdfj l ~u fl e + u e ; in this case further acceptance cuts have been imposed 
on the charged leptons: E\ > 20 GeV, px, > 10 GeV , \rji\ < 3. We see that the separation 
between the two Higgs mass hypotheses persists also in this more realistic setting. 

As an example of the different kinematical distributions in the two cases we show in 
Fig. [7] the absolute value of the difference between the pseudorapidities of the two Ws. 
It is interesting that the shape of the kinematical distributions, which are less sensitive to 
pdf uncertainties than their absolute normalization, behave differently when a light Higgs 
is present in the spectrum then when its mass is very large. For more details we refer to 



261 ]. We conclude that while it does not appear to be possible to study the contribution 
of the scattering diagrams in isolation from the remaining ones, the full amplitude, in the 
region of large WW invariant masses at LHC energies, is sensitive to the presence of a light 
Higgs and therefore to the details of the mechanism of EWSB. 

VII. CONCLUSIONS 

We have critically examined the role of gauge invariance in W-fusion processes and the 
reliability of the EVBA in describing them in the Unitary, Feynman and Axial gauge. We 
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FIG. 7: Distribution of the pseudorapidity difference between the two W's for us — > dcW + W at 
the LHC for the Infinite Higgs mass case (solid curve) and M/ l =200 GeV (dashed curve). 



have shown that the WW scattering diagrams do not constitute the dominant contribution 
in any phase jspace region for the set of gauge fixing we have examined. The Axial gauge as 
proposed in [12[| results in less severe cancellations between the contribution of non scattering 
diagrams and the contribution of the WW scattering ones but typically the two sets have 
comparable magnitude. 

We have shown that EVBA results and their relationship to exact results depend quite 
sensitevely on the set of cuts which need to be applied in order to obtain a finite result. 
Therefore, it is extremely difficult to extract from the EVBA more than a very rough estimate 
of the actual behaviour of the Standard Model predictions for boson boson scattering. 

We conclude that while it seems impossible to isolate the contribution of the scattering 
diagrams, the mechanism of EWSB can be investigated, using the full amplitude, by a careful 
analysis of the region of large WW invariant masses at the LHC. 
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APPENDIX A: AXIAL GAUGE 



We report here for convenience the main formulae describing the Axial gauge formulation 



of 



121 ] to which we refer for a more detailed discussion. 
Parametrizing the Higgs field as 

/ 

<j>(x) 



iw(x) 



i y \[v + h(x) + iz{x)\ 
the propagators for the five dimensional vector field (W^, w) is 



^ u = ( A1 ) 

y lvl w 

= - r + 551^ _ _5>v ( A2) 

g • n (g • nj" 1 

JV»(,) = -Ut w <-™" + v*' (A3) 
(g • n)^ 

N sv (q) = N us (q)* = N vs (-q) (A4) 
M 2 n 2 

*-(,) . 1 - ^ (A5) 

and the index s indicates the scalar component. 
The polarization vectors satisfy 

N IJ (q)= £ e'(g,A)e%A)* (A6) 

X=1,2,L 

For A = 1, 2 they describe transverse polarization: 

e s (q, A) = q^(q, A) = n^(q, A) = e M (g, A)e"(g, A')* = $aa' (A7) 
while for A = L they describe longitudinal polarization: 



where 



e s (q, L) = -iy/l- M^tf/iq ■ nf (A8) 

= -[M w /q-n}n^+[M^n 2 /{q-n) 2 }q^ 
yJl-M&n 2 /(q-n) 2 



The propagators and polarizations for the (Z^,w) field can be obtained from Eqs. (IAlHA9p 
with the substitution — > M/. 

The remaining Feynman rules are identical to the ones in gauges. In our calculation 
we have used £ = 1. 
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